UNCLASSIFIED 


r 


(^efvtaduced 
luf,  tt» 


ARMED  SERVICES  TECHNICAL  INFORMAIIQN  AGENCY 
ARLINGTON  HALL  STATION 
ARLINGTON  12,  VIRGINIA 


UNCLASSIFIED 


NOTICE:  When  govenment  or  other  drawings,  speci¬ 
fications  or  other  data  are  used  for  any  puxpose 
other  than  In  connection  %rlth  a  definitely  related 
govemoent  procurement  operation,  the  U.  S. 
Government  thereby  Incurs  no  responsibility,  nor  any 
obligation  idiatsoever;  and  the  fact  that  the  Qovem- 
ment  nay  have  fomulated,  furnished,  or  In  any  way 
supplied  the  said  drawings,  specifications,  or  other 
data  Is  not  to  be  regarded  by  implication  or  other¬ 
wise  as  In  any  manner  licensing  the  holder  or  any 
other  person  or  corporation,  or  conveying  any  rl^ts 
or  pemlsslon  to  manufacture,  use  or  sell  any 
patented  Invention  that  may  In  any  way  be  related 
thereto. 


S  D  A 


1 


Inf7,.  r».  Formnlo  4  Ing.  A.  minr-llnl 


.p.A.,  Via  Tlbartlna  Kir.  IP.'iOo,  ROMh,  ITALY 


TECHNICAL  ^‘'.CIF.NTIFTC  >  NOTE  NO.  ? 
CONTRACT  NO.  AF  6l  '05?^  -  “''*5 

r.ey  ^ 


Thp  TPSpar  (’h  reported  In  this  dorument  has  bern  sponsored 
by  the  Air  Force  Cambrldpic  Research  Center  of  the  Air  Re¬ 
search  and  Pevelopment  r'ammand,  United  States  Air  Force, 
thro'iph  Its  European  of  lice. 


(oqho 


di 


S.  p.  A. 


Abstract 

The  propagation  characteristics  of  a  uniform  magneto- 
ionic  duct  of  circular  cross-section  are  determined  for  fre¬ 
quencies  higher  than- the  cyclotron  frequency-  The  ratios  (plas¬ 
ma  wavelength)  /  (free  space  wavelength)  and  (power  flowing  In 
the  plasma)  /  (power  flowing  outside)  are  evaluated  and  discussed 
as  a  function  of  the  diameter  /  wavelength  ratio  and  of  the 
plasma  permettlvlty  for  the  propagation  of  circularly  symme¬ 
trical  modes.' 

The  Brlllouln  diagrams  for  frequencies  higher  and  lower 
than  the  cyclotron  frequency  are  also  derived. 
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Introduction 

_ * 

> 

In  Technical  Note  N.1  we  have  Investigated  electro¬ 
magnetic  wave  propagation  In  a  magneto-lonlc  duct,  namely  In 
a  cylindrical  plasma  of  circular  cross-section  In  a  static 
axial  magnetic  field,  assuming  signal  frequencies  lower  than 
the  cyclotron  frequency. 

In  this  Note  the  same  general  analysis  Is  extended  to 
the  case  of  signal  frequencies  higher  than  the  cyclotron  fre-  . 
quency.  Here  too  the  discussion  Is  based  on  the  general  sha- 
raaterlatla  equation  and  restricted  to  a  uniform,  non-atten¬ 
uating  plasma  and  to  circularly  symatetrlcal  modes:  the  permet- 
tlvlty  tensor  Is  eftosen  as  a  constant  parameter  for  each  solu¬ 
tion.  From  these  and  the  previous  results  we  have  derived  the 
Brlllouln  dlagraais  relating  the  guided  and  the  free  space  phase 
constants  for  a  given  set  of  plasma  frequency  and  cyclotron 
frequency  values.  Tltese  diagrams  are  similar  to  those  derived 
with  the  so-called  quasl-statlc  approximation,  except  for  a 
few  special  features,  which  will  be  described. 


$ 
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1)  Sunmar^r  of  th#  previous  results  t 

The  dielectric  constant  of  a  uniform  plasma.  In  the 
absence  of  agitational  or  drift  motion,  placed  In  a  steady- 
state  magnetic  field  B  directed  along  the  positive  r(l,l  -  ?) 
axis.  Is  a  tensor  //<f//  defined  by  the  formula: 


where 


-P.-  = 


<*  /'•**  /  *•  •: —  •  — r~  * 


*  fO 


A/- 


‘^1.' 


»e 

«  & 


n  -  electron  density 
m  >  electron  mass 
e  «  electron  charge 

From  the  equatlonsCi ) ,  one  can  easily  find  the  follow¬ 
ing  relation  among  ^  and  fg 

e\  •  (t) 


If  the  plasir.a  Is  uniform  and  In  steady  state  conditions 
and  we  take  a  system  of  cylindrical  orthogonal  coordinates 
f,fi^  where  the  z  axis  Is  parallel  to  B  and  with  the  same 

f 

orientation,  the  Maxwell  equations  can  be  solved  and  the  field 
equations  for  the  circularly  symmetrical  modes  written  : 
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Inside  the  plasma: 


*  A«  L,  i,  ^ 

f^*  [fi,Aj.(x.f/ie)tfit  ^  J,r^^f/^)J  (3) 

H.,.iES±  [ft.lJ,(^.r/^)-f  ^  rx,//£)J 

,  .  j  «<  i  [fl,  ^  *  fl.  i  J,  (^.  f/t)  J 

•\9 


outside  the  plasma: 


where ; 


B,  s  &  Ko  {x,  ) 


./i^  oc  6  -L  K.rx.f/fi) 
Ao  '*‘® 


t  y7  - 

H,  > 
hi  t 
Hi 


.  C  JL  K,fX(,jo//e) 

ho  *-» 

.  2^  c  -L  K/^of/e j 

Xo 

*  -•’^  i;  *:,<•*•/’/«) 

=  j4-  •  c  ■  ) 

c  ^ 


R«  S  =  radius  of  the  plasma  cylinder 

j/ . ^  .  free  space  to  guided  wavelenght  ratio  (  >  l) 

Z,"V|^-  free  space  characteristic  Impedance 


foglio  ^  di  '*'9 


S.  p.  A. 


Morieover ; 


and 


AT,,.  -  iii 


(5) 


y«  - 


“TTiCo"  ■"  a£,  [  Y  rv‘-</  - 


'*0 


0<y.. 


If  the  boundary  conditions  are  Imposed,  the  following 
equations  result: 


Ji  •*•  W, 


X*  f  ) 

a.  7,u)  *3.^  JotxQ 
K#  r  X*; 

' 

j.r^o 


»•  1 _ 

h  K,  (^*) 

T,  ^•(M 
».  T/.V  K, 


Ko(>‘^'> 


The  last  two  equations  are  consistent  only  If 


F  (x^  ,)]  ,  =  r  F(k,  ^  (7) 
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where 


For  each  assumed  set  oIV,t,  and  values,  the  variable  A'o  Is 
the  only  unknown  In  (7'S  and  It  can  be  determined. 

The  value  which  corresponds  to  the  assumed  u  , 
and  set  la  then  computed  from  the  equation: 

i  x,^TrV(K^-i) 

C 

When  S^<o  and  oc  Is  smaller  than 

oC.  =  (s.-fi)  3 

the  quantities  are  complex. 

In  this  case  we  write: 

»^p(  tije) 

«  eS^f/rix-p,  9) 

so  that  equation  (7^  can  be  written  as  an  equation  of  only 
real  quantities  as  follows: 
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t  %  H|  ( *’») 

1^1 


_l_  ^  ^ 

T  KoNI-  "■-' 


AAM  i 


[(^-ff)HV.-Y,)J 


J 


The  parameter  can  thus  be  computed  solving  this 

last  equation.  ^ 

It  is  convenient  to  divide  the  entire  range  of 
and  £i  values,  which  may  exist  according  to  eqs .  fl),  into 
four  regions: 


1)  I  ^  *0O 

2)  0<«,<l 

»  • 

3)  .ft?  <  S';  \  0 

0<  fi<  I 

4)  .»<£,<  i> 


The  first  region  was  examined  in  the  Technical  Scien¬ 
tific  Note  N.  1.  All  the  other  tjirec  regions  satisfy  the  con¬ 
dition  to  >  and  are  examined  here. 


joglio  ^  di  ^9 


S.  p.  A. 

2)  Propagation  in  the  o<f, <)  ^  region 

dn  this  /'eglon  the  angular  frequency  ^8  higher  than 


4- 

.  A  simple  discussion  of  eqs.  (6)  and  (8)  shows  that 
r  ^  i  ‘ 

3,, I  and  T  arc  always  negative  for  .  Since 

r  Is  then  always  negative  and  positive  (sec  flg.O*  r.o 

Solution  exists  In  this  region.  The  conclusion"  la  that 

Is  a  cut-ofi  frequency  and  no  propagation  is  possible  forw^^^f 


"5^  Propagation  in  the  •oo<£,<o,  o<  <  i  region 

f-i 

Here 't'  and  are  always  negative  (see  fig.  2,  3)  and 

b,  positive.  The  curves  of  fig.  1  are  stlli  valid  In  this 

case ‘and  show  that  an  Infinite  s^et  of  solutions  exists. In  the 

/ 

frequency  range  this  region  corresponds  to  the  conditions 
^  <  w  <.  vJ  ^  If  and  If 

When  o^^approaches  unity,  y*  tends  to  zero  and  di/X„ 
to  one  of  the  values  (el/Ao)^ 

J<^) 


where  Is  the  «  root  of  J 


When  approaches  Infinity,  'T'  tends  t.  zero 
and  the  solution  tends  to  one  of  the  poles  of  P^,  namely 

to  Too ,  where  Slor-f^ie,  and  Ac* ,f*)  Is 

A.  W 


the  m^*^  root  of  the  equation 


J,  ) 


Jo(^o»X'o) 


K|  (Xe] 
Ko  (^o) 
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Consequently  oL^Xo  approaches  zero. 

The  behaviour  of  Is  shown  In  fig.  4  for  the  special 
case  ■  -2,  «  0.5,  while  In  fig.  5  and  6  we  have  plotted 

^  ^and  versus  i/X,  (  P^^nd  being  respectively, 

the  power  flowing  Inside  and  outside  the  plasma). 

The  very  unusual  characteristics  of  this  type  of  pro¬ 
pagation  are  summarized  below: 

1)  For  each  mode  propagation  can  take  place  orily  for 
values  of  the  ^/Ao  ratio  between  zero  and  a  maxlrnum 

which  Is  characteristic  of  that  mode. 

2)  For  a  given  mode,  wher^  d/Ap  is  between  f<A/Ao')^ 

and  propagation  can  take  place  with  two  different 

Aj/A.0  ratios,  but  their  group-velocities  ere  opposite  In 
sign. 

51  The  axial  flows  of  electromagnetic  power  In  the  air 
and  In  the  plasma  have  generally  opposite  dlrsstions.  If  Pp^f^ 
we  have  an  "average”  forward  wave.  If  Pp  <  Pg  an  "average" 
backward  wave.  The  points  labeled  "A"  In  fig.  5,  where  the  group 
velocity  Is  zero,  correspond  to  the  points  "A"  In  fig.  6  where 
Pf- -  Pe»  ,  so  that  although  some  power  Is  really  flcwlng  sep¬ 
arately  In  the  air  and  In  the  plasma,  the  total  average  flow 
Is  zero. 

4^  The  highest  ratio  Pp /  is  obtained  at  ol/Ao=OJ 

at  this  point  the  m*^  mode  power  ratio  Is: 


^  K.qT)  _  K,V^) 

4)  Propagation  In  the  <0  ,  region 

Here  In  the  <<;<oC  <oo  range  Is  srs  negative 

flg.7>  and  'P  Is  positive  and  smallsr  than  unity^flg.  B)  ; 
the  behaviour  of  P,  ..  for  the  general  case  Is  shown 

I  f  tr  * 

Ir.  fig.  o. 
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The  curves  indicate  that  in  tnla  region  only  oue  Xo  solution 
may  exist. 

m  / 

4 

When  oL  approacxies  infinity,  y  tends  to  zero  and  the  X© 
solution  of  tne  dispersion  equation  tends  to 

Tne  parameter  X*  increases  wnen  o<  decreases  and  approacnes  in¬ 
finity  at  ^  ,  tne  root  of  tne  equation; 


when 

ir,i 

iTl’/i 

when  <  <  ^ 

iri 

Wnen^/oi-  t«e  X«  solution  crosses  continuously  the  rn- 

abscissa  and  tends  to  infinity  in  tne  complex  region  . 

Curves  of  are  plotted  in  figs.  V)  to  tor  tne  spe- 

•  ^ 

clal  case  -  -2,  -  -  1.5,  while  the  X,  versus  U  curve  is 

given  in  fig.  1^. 

When  tne  curves  snow  clearly  tnat  no  solution 

exists  at o^jtnen  we  do  not  expect  any  solution  also  at  finite 
values.  Ine  propagation  region  is  thus  bounded  by  tne  condi¬ 
tion  -5, 1  and  becomes*. 

f  ,  <£.■£-  I 

I  £, 

In  tne  frequency  domain  this  corresponds  to  tne  range: 

1*/^  .c  tv  < 

Wnit  tne  usual  formulas  tne  \^/L  and  versus  d/lo  curves 
cem  be  calculated;  tney  are  plotted  in  figs.  H  and  1';,  for  tne 

I 

case:  *,  =  -2;  =  -1.'-'. 

We  snail  finally  note  tnat,  as  in  tne  previous  case,  tne 

power  outside  the  plasma  flows  in  a  direction  opposite 
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to  that  of  power  inside  so  that,  being  generally  fp  « fp  , 

« 

the  power  flow  in  the  plasma  is  backward  although  the  total 
power  is  flowing  forward. 

BRILLOUIN  DIAGRAIC 


Keeping  <f,  and  i-'j  as  constant  pa-^araeters  in  the  Max¬ 
well  equations,  it  has  been  possible  to  solve  the  dispersion 
equation  with  a  straightforward  procedure.  Brillouin  dia¬ 
grams  (  versus  )  are,  instead,  more  difficult  to 

derive  because  e,  and  are  functions  of  frequency. 

In  order  to  obtain  from  our  computed  data  these  dia¬ 
grams,  we  have  to  solve  the  following  system: 

X./Ao  =  -T  ( •^/Ao  '  ,h) 

t,  -  (  /  .  4^-.  ) 

'  (1.4.) 

lc‘  I  ' 


(10) 


Let  us  keep  constant  the  rati  o  oi --  and  choose  ail  tne  pair^ 
of  and  values  which  satisfy  the  equation: 


(11) 


0.0 -fi) 


( see  fig.1 6 ) 


The  first  two  equations  of  the  system  (10)  allows  us  to  plot 
in  the  versus  plane  a  set  of  curves,  i'j 

being  a  constant  for  each  curve.  Let  us  find  the  intercepts 
of  these  curves; 


(liO 


K  i  =  ^  {  'h  d.,  0^ , 


With  the  straight  lines; 
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htnich  represents  tne  tnlrd  equation  of  tne  above  system 
(10)«  Kpd  being  tne  plasma  density  parameter  - 

For  given  values  of  the  frequency  Independent  para¬ 
meters  a  and  Kpd«  these  intercepts  at  different  values 
plotted  In  tne  Kd  versus  plane  provide  the  Drlllouln 
diagrams: 

(14)  ^  =  FC  KM  ,  d.  ,  KpM.) 

The  computation  nave  been  performed  for  two  a  values, 
tne  first  being  chosen  In  tne a>t region  (as  2.25)»  tne  second 
in  thea.<i  regloji  (a  ■  0.5).  la  figs.  1?  to  20  we  nave  plotted 
tne  curves  given* by  tne  equation  (12)  for  the  first  X*  solutions 
and  respectively  for  tne  conditions  I)  a  =»  2.25, 

II)  a  s  2.25,  III)  a  =  0.5,  IV)  a  -  0.5  -7^0  - 

Tne  Brlllouln  diagrams,  wnicn  result  from  tnese  and 
similar  curves  for  tne  second  X.  solutions,  are  plotted  in 
figs.  21  to  24. For  eacn  case  three  Kpd  values  nave  been  consi¬ 
dered,  namely:  I.5,  5.5,  10. 

It  is  interesting  to  compare  these  curves  with  the  corre¬ 
sponding  curves  obtained  by  Could  and  Trlvelplece^^^and  ^>y  Smul- 
lln  and  Cnomey^^^uslng  tne  so-called  quasi-static  approximation. 
These  authors  i^ave  studied  tne  more  general  case  of  a  circular 
waveguide  partially  filled  with  plasma;  from  tneir  results  tne 
curves  to  be  used  in  our  case  are  easily  derived  by  setting  tne 
waveguide  diameter  equal  to  infinity. 

No  significant  difference  is  found  in  tne  benavlour  of 
tne  Brlllouln  diagrams  in  the region  (Kd<  K^). 

It  this  region  all  tiie  Kd  versus  Sd  curves  start  from  tne 
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/  origin  and  approacn  asymptotically  tne  value 

or  ,  whichever  is  smaller,  except  the  first  solution 
mode  which,  when  «  approaches  the  asymptotic  value 

A  behaviour,  slightly  more  complicated  than  predicted  from 
the  quasi-static  approximation,  i  s  found  in- the  .>  u>p 
region.  Here  the  curves  start  1 rom  different  points  of  the 
Kd'jid  line  between  Kp«l~and  and  approach  the 

asymptotic  value  ^^4  or  K^ct  ,  whichever  is  larger.  The  slopt.  ^ 
of  these  curves  may  bo  positive,  negitive  or  partially  po-  ^ 
sitivc  and  partially  negative,  so  that  tie  waves  can  be 

forward  or  backward,  it  Is  worth  to  not*'  that  ir,  th'-  nobs' 
static  appr  ox  Ima  t  Ion  nl’.  tnes'  wnvi»s  are  always  h-r^war  rt  . 

The  following  details  can  be  added  to  the  discussed 
fent’ires  in  the  region .  Tlie  irillouin  diagrams  start¬ 

ing  points  on  the  Kd  -/Jd  Un'^  are  glv>*n  from  th«'  simple 

relation:  , - - 

K  4  .  K a 

I  *■  c 

where  :  c  =  2 

Anofner  simple  relation  can  be  found,  which  allows 

U3  to  distinguish  between  forward  and  backward  wavfs  in 

• 

the  ^  region,  l^et  uc  observe  th»  behaviour  of 

the  intercepts  which  provide  the  orillouin  diagrams,  when 
it  increases  from  zero  to  unity  and  Kd  varies  consequently 
from  to  infinity.  indicates  clearly  tiiut,  when 

the  intercept  points  fall  in  toe  positive  S,  region,  t.iere 
io  only  one  intercept  for  each  value.  In  this  case  the 
Kd  versus  ^4.  curve  has  cverj-vhere  a  por.itivc  slope, 
tite  limit'  ^d-^oo  will  be  attained  wuen  Kd-  since  It  nas 

to  be  considered  as  t!ie  intercept  with  the  curve, 

which  is  the.  only  curve  (id)  witli  an  horizontal  asymptot 
(Kd  -  a  /(■  ft.  - 1 )  ) 
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A  s.traitforward  discussion  t;ien  shows  that  all  the  inter- 
cept  points  of  a  Briliouirt  curve  shall  fall  in  the  posi¬ 
tive  region  as  far  as  is  larger  than  this  asym- 

potic  value  .  This  is  thus  the  condition  for 

forward  waves  and  cc*  b*  rewritten  in  tne  form: 

When  this  condition  is  not  satisfied,  the  waves  may  be 
backward  or  start  at  low  /icL  values  as  forward  and  change 
to  backward  a \  large  values. 
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